Two-dimensional electron systems beyond the diffusive regime 



O 
(N 



P. Markos 
Department of Physics FEI, Slovak University of Technology, 81219 Bratislava, Slovakia 

Transport properties of disordered electron system can be characterized by the conductance, Lya- 
punov exponent, or level spacing. Two additional parameters, K\\ and 7 were introduced recently 
which measure the non-homogeneity of the spatial distribution of the electron inside the sample. 
For the orthogonal, unitary and symplectic two dimensional disordered models, we investigate nu- 
merically the system size dependence of these parameters in the diffusive and localized regime. 
Obtained size and disorder dependence of Kn and 7 is in agreement with with single parameter 
transport theory. In the localized regime, 7 — > independently on the physical symmetry of the 
model. In the diffusive regime, 7 equals to the symmetry parameter /3. For the symplectic model we 
analyze the size dependence of 7 in the critical region of the metal-insulator transition and found 
the non-universal critical value y c - 



PACS numbers: 73.23.-b, 71.30.+h, 72.10.-d 
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I. INTRODUCTION 

Transport of electrons through disordered struc- 
tures offers a broad variety of interesting universal 
phenomena*^ With increase of the strength of the dis- 
order the character of the transport changes from the 
ballistic to diffusive up to the insulating, where all elec- 
trons are localized^. 

In the limit of weak disorder (diffusive regime) the 
transport can be studied analytically using, for in- 
stance, the Dorokhov Mello Pereyra Kumar (DMPK) 
equation^ the Green's function analysis^ or random ma- 
trix theoryj&i The existence of the metal-insulator tran- 
sition in two and three dimensional models^ is a strong 
motivation to construct an analytical theory of the trans- 
port beyond the diffusive rcgim o 10 ' 11 . Also, numerical 
data for the localized regime 1 ^— show that, contrary to 
theoretical expectation, the distribution of the logarithm 
of the conductance is never Gaussian for disordered sys- 
tems in higher dimension. Therefore, a general transport 
theory must explain how the dimension of the system and 
physical symmetry of the model^ influence the ability of 
electron to move through the sample. 

The most elaborated analytical description of the 
transport in strongly disordered structures is based on 
the generalized DMPK equation (GDMPKE) £L The the- 
ory takes into account that the spatial distribution of 
electrons in the regime of localization is not homoge- 
neous. The last was confirmed by numerical simulations 
in Ref . 14 i 18 In GDMPK, the non-homogeneity of electron 
distribution is measured by a large number of parameters 
K a b (defined later); however, only two of them, Kn and 
7 = 2A'i2/i\n are decisive for the transport^ 

The GDMPKE is not exactly solvable, but approxi- 
mate analytical solution for 3D disordered systems^— 
agrees very well with numerical data. Numerical solu- 
tion of GDMPKE^ 2 - confirmed that it correctly describes 
disordered orthogonal systems and that parameters K a b 
depend on the dimension of the system. 

Detailed numerical analysis of parameters An and 7 
in three dimensional model was performed in2£. The aim 



of this paper is to investigate how these parameters de- 
pend on the physical symmetry in two dimensional (2D) 
models. We present numerical data for the parameters 
An and 7 for the orthogonal model (O), unitary (U) 
and two symplectic (S) 23 i 24 models in diffusive and in- 
sulating regime. For the S models, we also study the 
behavior of both parameters in the critical regime of the 
metal-insulator transition. 



II. GENERALIZED DMPK EQUATION 

Consider a disordered system of the length L z con- 
nected to two semi-infinite ideal leads with N open chan- 
nels. Transmission parameters are given by the transfer 
matrix, which can be written in general form as^ 



T = 



u 
v! 



VTTx Vx 
Va vT+a 



v 
v' 



(1) 



In Eq. (|TJ|, u, v are N x N matrices, and A is a diagonal 
matrix, with positive elements A Q ,a = 1,2, ...N. In sys- 
tems with time reversal symmetry, matrices u' and v' can 
be represented in terms of u and v^ For the orthogonal 
system, vl = u* and v' = v* . For the symplectic symme- 
try, the scattering depends on the spin of the electron; 
the elements of matrices u and v are 2x2 matrices which 
fulfill the symmetry relation o 6 ' 25 
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Statistical variables u, v and A contain entire information 
about the transport. In the weak disorder limit^ the con- 
ductance g (in units of 2e 2 /h) is completely determined 
by eigenvalues A a £i2£ 
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N N 

*-i 1 + A a ' = ^ 



1 



^ cosh x a /2 



(3) 



In the last equation, we used the parametrization A Q 
(cosha; a — l)/2. 



The probability distribution of As can be found as a 
solution of the DMPK equation^ The generalization of 
the DMPK for the orthogonal symmetry class, was done 
by Muttalib and Klauder— who introduced new param- 
eters, K a b which characterize the spatial distribution of 
the electron in the disordered sample. The generalized 
DMPK equation rcadsil 
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where £ is the mean free path, and 
N 

a<b 
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This equation can be simplified when all K aa are ap- 
proximated by K\\ and ~f a b ~ 7 for all a, b (a ^ b). This 
approximation was confirmed by numerical wor k 20 ' 22 . 

Although the conductance is still given by Eq. ((3j), it 
becomes implicitly a function of the spatial distribution 
of the electron. 



III. MODELS 

In numerical work, disordered sample is represented by 
two dimensional (2D) square disordered lattice of the size 
Lx L. The orthogonal 2D model with on-site disorder is 
defined by the Hamiltonian 

H= W^xzCxzclzCxz 

+ 1/ lE IZ C i+a/^+4 C x+^ (6) 

' '|| Z-txz C x,z+a c xz + C xz C x ^ z +a 

Here, a is the lattice spacing, e xz are random energies 
from the box distribution, \e xz \ < 1/2, W measures the 
strength of the disorder and Vj| = 1 defines the energy 
scale. To avoid closed channels in leads, we use V±/Vu = 
t < li^i In what follows we consider t = 0.9, the energy 
of the electron E = 0.01. With a = 1, wc identify the 
number of channels 



N = L. 



(7) 



It is generally accepted^ that only localized regime 
exists in the model when the size of the system L — > oo 
(the critical disorder W c = 0). Nevertheless, diffusive 
transport is observable for sufficiently weak disorder and 
small sample size»^ 

The second model of interest is the symplectic model 
with spin dependent hopping. Here, the hopping of elec- 
tron from one site to the neighboring one can be accom- 
panied by the change of the sign of the spin and Vji , V± 
become 2x2 matrices. In numerical simulations, we 
study the Ando model with hopping hopping terms 



V±=t 



Vi -v 2 



V, = 



Vi 



-iVo 



-iV 2 Vi 



(8) 



The spin-orbit coupling is characterized by the parameter 
S = Vi and V^ + V^ 2 = 1. In this paper, S = 0.5. We also 
study the Evangelou-Ziman (EZ) models which uses the 
random hopping matrices V: with help of three indepen- 
dent random variables, t x ,t y ,t z , distributed uniformly in 
interval (— /x/2,/x/2) 

V -V - f f 1 + itZ ~t v +it x \ (q , 

Vx-V xz , x+az -ti tv _. tx 1 _ ifZ , (9) 



and 



*\\ *xzxz-\-a 



1 + it z -t v + it x 
t y -it x 1 - it z 



(10) 



and consider (i = I. 

Both Ando and EZ model exhibit the metal-insulator 
transition when the disorder W reaches the critical value 
W c i 23 i 24 Owing to the anisotropy of our models, the 
critical disorder differs from that obtained in previous 
works. 23 ' 28 We found W c sa 5.525 for the Ando model 
and W c = 6.375 for the EZ model. 

The 2D model with external magnetic field B can be 
obtained by including the Peierls hopping term V± = 
texpixa, a = (e/K)Ba 2 into the Hamiltonian ([6]). 



IV. THE MATRIX K 

The matrix K a f, is defined in terms of higher moments 
of the matrices v. 

Kab = (kab) (11) 

Here, (...) represents an ensemble average. 

For the orthogonal system, the matrix k® b is defined 

„„17 



k ab = ^\v a a\ 2 \v a b\ 2 '■ 



In the diffusive regime,- 
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For the unitary models, 



1 + Sab 

L + l 



fc Q u b = El w -l 2 Ko 



and 



K u — — 

1^ 



L a b 



(12) 



(13) 



(14) 



(15) 



For the systems with symplectic symmetry the matrix 
° b is given. 



fc ab = X! V ia v * a bV a bV a 
a=l 



(16) 



In this equation, the 2x2 matrices v^ , v* and v are 
defined in terms of the matrix v: 



and 
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In the diffusive regime, if| 6 is degenerated diagonal 
matrbs^E with diagonal elements 



A a6 — 



2 - <?>ab 

2L-1' 



(19) 



Our numerical results discussed in Sect. [VI confirm that 
the same holds for any disorder strength. 
From Eqs. (|13I15I19[) it follows that 



E^°=E< 



J2 K * = 1 - 



(20) 



We use these relations to test the numerical accuracy of 
our results. 
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FIG. 1: (Color online) K\i and K12 as a function of the 
system size for square lattice L x L for the symplectic (S) 
and orthogonal (O) s and unitary (U) systems. The disorder 
W = 2. Solid lines are linear fits with slopes given in the 
legend. 



Obtained data for k a \, were also used for the calculation 
of probability distributions. 

As noted in Section ITVl K„ h are 2x2 matrices. Nu- 



merical data confirm that, with the relative accuracy of 
1CP 3 , these matrices remain diagonal degenerate for each 
value of the disorder and all size of the system. 



We consider square samples of the size L x L attached 
to two semi-infinite ideal leads. The size L increases 
from L = 10 to L = 256 (S model) up to L = 600 
(O model). For each value of L and W, we analyze 
the statistical ensemble of typically iV s t a t ~ 10 4 sam- 
ples (iVfjtat ~ 1000 — 4000 for the largest system size). 
In numerical calculation, the sample and leads are rep- 
resented by the 2 A x 2 A transfer matrix M and Mo, 
respectively^ Following^ 2 -, the conductance is given as 
a trace of matrices L + MR + and L~MR~, where i?" 1 
(L+~) arc N x 2N (2N x N) matrices composed of right 
and left eigenvalues of Mo, respectively. The upper index 
+ (~) indicate the direction of the propagation through 
the sample. Comparing with Eq. ([T]) we find 



L+MR+ = v(l + A) 



-v 



and 



L~MR- =v' r (l + \y 



(21) 



(22) 



A. Diffusive regime 

We first verify the prediction of the DMPK equation 
for the diffusive regime. In Fig. [T]we show the L depen- 
dence of parameters K\\ and K12 for the orthogonal and 
symplectic system with disorder W — 2. The system is 
in the diffusive regime (the conductance g varies between 
4.9 and 5.03 for the orthogonal model, and increases from 
7 to 11 for the S model). Linear fits shown by solid lines 
confirm that both K\\ and K\i ~ \/ L and 7 equals to the 
symmetry parameter (3 in the diffusive regime. The spa- 
tial distribution of electrons is homogeneous and no addi- 
tional parameter must be introduced into the model. The 
transport is universal, the only model parameter in the 
DMPK is the ratio L z /£ of the system length to the mean 
free path. Although the DMPK was derived only for the 
quasi-one dimensional systems, our data show that rela- 
tions (|13ll5ll9p are valid also for the square samples. 



Thus, eigenvalues A a can be obtained numerically by di- 
agonalizing of the matrices L + MR + and L~MR~. Ma- 
trices v and v' consist of corresponding eigenvectors. De- 
tails of numerical method are given in Rcf^ Mean val- 
ues, Kn and A'12 were calculated as an average over the 
statistical ensemble 



K ab = 



1 



w., 



As, 



Y.K 



(23) 



B. Insulating regime 

In the limit of strong disorder, we expect that K aa 
depend on the index a and K aa ~ 0(1). Contrary, off- 
diagonal elements K a b, a ^ b, should decrease to zero, 
K a b ~ 1/i (a / !») so that 7 ~ 1/L decreases to zero 
when the size of the system increases^ 

Figure [2] shows the L dependence of A° and K^ for 
orthogonal systems with various strength of the disorder. 
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FIG. 2: (Color online) The size dependence of Kn (top) 
and if 12 (bottom) for the 2D orthogonal model for various 
strength of the disorder (given in the legend of bottom panel). 
A"ii converges to the non-zero limit when L — > oo for each dis- 
order strength. This limiting value, however, is too small to 
be observable numerically for weak disorder within the con- 
sidered size of the system. A12 decreases to zero for any value 
of the disorder W. 



Similarly to the 3D orthogonal model discussed ir>2£, both 
K^[ and K® 2 are linear functions of 1/(L + 1), Since no 
metallic regime exists for the non-zero disorder, we ex- 
pect that Kyy converges to the nonzero value in the limit 
of L — > 00 for all values of W . 
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The limiting value K ^ can be easily calculated numer- 
ically for strong disorder. This is more difficult for weak 
disorder (W < 4), since K® x becomes smaller than the 
inverse of the accessible sample size. 

Similar data (not shown) were obtained also for the 
symplectic models. 



C. Critical regime (symplectic models) 

Critical regime exists only for the S systems. In the 
critical regime, W = W c we found that both Kf x and 



FIG. 3: (Color online) The size dependence Kn and K\2 at 
the critical point for the symplectic (Ando and EZ model) 
models. Solid and dashed lines are power fits K\ a oc L~ a 
(a = 1, 2) with the exponent a « 1.003 - 1.005. 
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FIG. 4: (Color online) The size dependence of the parameter 
7 = 2Kf 2 /Kfi for the Ando model. In the metallic regime, 
(W < W c = 5.525), 7 slowly increases when L increases and 
converges to its isotropic value 7 = 4. In the localized regime 
(W > W c = 5.525), 7 decreases to zero when the size of the 
system increases. At the critical point, 7 does not depend 
on L. The L-independent critical value of 7 s m 2.60 (dotted 



(.") line) is plotted by dot line. 



Kf 2 decreases at the critical point to zero 
K n (W = W C ), Kf 2 (W = W c )oc 



L 



(25) 



(Fig. [3]), so that 7 s reaches a critical value, 7° = 
2Kf 2 /Ki 1 which does not depend on the size of the sys- 
tem 



7 C = const. 



(26) 



As shown in Fig. [31 the critical value 7^ is not universal 
but depend on the model. We obtain 7 s = 2.601 for the 
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FIG. 5: (Color online) The size dependence of the parameter 
7 = 2K° 2 /KrL for the 2D orthogonal model. We expect 
that 7 decreases to zero for any value of the disorder W. 
However, the size of the system L is not sufficient to see this 
decrease when the disorder is small. Instead, we observe the 
diffusive limit 7 ~ 1. Only for sufficiently strong disorder, 
7 decreases to zero. In contrast to the symplectic model 
(Fig. 3]) there is no indication for the existence of the critical 
point where 7 converge to the L-independent limit. Inset 
shows the disorder dependence of 7 U for the unitary ensembles 
(L = 128, a = 1/8). 



Ando model, and 1.795 for the Evangelou-Ziman model. 

Figure U shows that the length and disorder depen- 
dence of parameter 7 s can be, at least in principle, used 
for the estimation of critical parameters in the same 
way as mean conductance of the smallest Lyapunov ex- 
ponent. For very weak disorder, we find that 7 s only 
weakly depends on the size of the system and increases 
to the metallic limit 7 = 4 when L increases to infin- 
ity, indicating that the system is in the metallic regime. 
For stronger disorder, 7 oc 1/L decreases to zero when 
the size of the system increases, in agreement with the 
prediction of the theory.— We found the critical regime 
between these to limits, where 7 converges to the size- 
independent constant 7 C ~ 2.60 (obtained already in Fig. 
when W = W c . 

For comparison, we show in Fig. [5] data for 7 calcu- 
lated for the 2D orthogonal model. Wc found no critical 
regime. Although 7 ss 1 for weak disorder, we expect 
that this is the finite size effect, and 7 will decrease to 
zero for each disorder strength when L increases £ 



D. The universality 

With two new parameters Ku and 7, we must ver- 
ify if the transport properties of the system are still 
maintained by only a single parameter^ In the metal- 
lic regime, the answer is trivial since the entire matrix 
K reduces to model-independent numbers given by Eqs. 





FIG. 6: (Color online) Ku as a function of the localization 
length £ (estimated as £ = 2L/xi) for the 2D orthogonal 
model (top) and the 2D symplectic model (bottom). The 
strength of the disorder varies between W — 1 and W — 16 in 
both Figures. For strong disorder (2L/xi — > 0) data converge 
to the universal curve (solid line is a function 0.288 — 0.066 ln£ 
for the orthogonal model). Insets in both panels show that 7 
is an unambiguous function of xi. 



(|13|15|19|) . The universality of the critical regime was 
shown in the previous section. Here, we concentrate on 
the localized regime, where we expect that Ku becomes 
an unambiguous function of the localization length^ 
The last can be estimated from the smallest parameter 



Xl, 



Xl 



(27) 



In Fig. [5] we plot Ku as a function of £ for the orthogo- 
nal and symplectic Ando models. Data confirm that the 
parameter Ku becomes an linear function of ln£ with 
increasing system size and converge to the system-size 
independent limit when L — > 00. 

Two inset of Fig. [51 show that the parameter 7 is 
an unambiguous function of xi in all three regimes. In 
the localized regime, when Xi ~ L, data confirm that 
7 ~ 1/L, consistent with prediction of the Muttalib's 
theory. 
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FIG. 7: (Color online) The probability distribution of fen 
for the 2D orthogonal model with various disorder strength 
(given in the figure) . The size of the system L = 200. Inset 
shows the distributions P(ln An). Data confirm that fen is a 
good statistical variable with a well defined mean value and 
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FIG. 9: (Color online) The probability distribution of fen for 
the 2D Ando model at the critical point W — W c = 5.525 
and in the metallic regime (W = 2, inset). Both the square 
samples 100 x 100 and Q1D samples 100 x 5000 are considered. 



regime. We demonstrate that the distributions for the 
square sample L x L with quasi-one dimensional systems 
are almost identical. 
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FIG. 8: (Color online) The probability distribution of fen in 
the strongly localized regime for the Ando model. The mean 
value Kn is given in the legend. 



E. Statistical properties of fen 



We have shown that K\\ — ¥ in the metallic regime 
but Kn ~ 0(1) in the insulator. Small values of Kn in- 
dicate that the mean conductance of the system is large. 
Contrary, large values of Kn correspond to systems with 
small mean conductance. This is in agreement with our 
expectation: small conductance means that the electron 
has problems to go through the sample. When it finally 
reaches the opposite side, its spatial distribution is not 
homogeneous any more. 

However, the correspondence large fen - small g holds 
only for mean values of these parameters. As shown in 
Fig. [TU1 the values of g and kn for a given sample are 
not correlated within a given statistical ensemble: small 
values g <C (.9) can be accompanied with any value of 
fen - cither small fen <C Kn or large kn ^> Kn- The 
absence of the correlation observed in both the metallic 
and in strongly localized regime, confirms that the sta- 
tistical fluctuations of kn do not affect the mean value 
of the conductance. 



In the previous analysis we dealt only with mean val- 
ues of A11 and K12. Since both kn and k\2 are sta- 
tistical variables, we must also to study their statistical 
properties. Figure [7] shows the probability distribution 
of parameter kn and In fen for the 2D orthogonal model. 
For each disorder, the mean value can be identified with 
the most probable value. In the localized limit, both Kn 
and var kn are of order of unity, and the distribution 
P(fcn) becomes size independent (Fig. [8J. 

In Fig. |H]we plot the probability distribution of kn for 
the symplectic Ando model in the critical and metallic 



VI. CONCLUSION 

The electron transport through disordered system is 
determined by spatial distribution of the electron inside 
the disordered sample, which can be measured by param- 
eters Kn and 7. Our aim in this paper was to investigate 
how these two parameters depend on the size of the sys- 
tem, strength of the disorder and physical symmetry of 
the model. We concentrated on 2D disordered systems. 
In order to better understand the role of the disorder, we 
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compare numerical data for the orthogonal and symplcc- 
tic physical symmetry. For completeness, we add also a 
few data for the unitary ensemble. 

In the diffusive regime, the size dependence of both 
parameters follow the analytical relations given by the 
theory of DMPK equation. In particular, 7 equals to 
the symmetry parameter /?. In the localized regime, K\\ 
converges to the size independent limit and 7 ~ 1/L. 

For the symplcctic models, which exhibit the metal- 
insulator transition, we analyze the size dependence of 
both parameters and we found that 7 possesses a crit- 
ical value 7 C when disorder W = W c . Also, we found 
no significant difference between the values of K for the 
two dimensional and quasi-one dimensional systems. No 
critical value was found for the orthogonal model. 

We also found that K\\ is an unambiguous function of 
the localization length £ and 7 is uniquely given by the 
parameter x\. Therefore, the use of these parameters 
does not contradict the single parameter scaling theory. 

Since the elements of matrices k are given by elements 
of statistical matrices v, they are also statistical vari- 
ables. Fortunately, analysis of their probability distribu- 
tions confirm that their mean values are good representa- 
tives of the statistical ensembles. We found no statistical 
correlations between the conductance and k\\. There- 
fore, we conclude that mean values, K\\ and K12, and, 
consequently, 7 = 2K12/ Kn, are physical parameters for 
the description of disordered systems. 



FIG. 10: 2D Ando model: fen as a function of the conduc- 
tance statistical ensemble of iVstat = 10 4 samples with disor- 
der W — 2 (top) and W = 9 (bottom) . The size of the system 
L = 100. 
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